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EVALUATION  OF  A  COMPARTMENTAL 
MODEL  FOR  PREDICTION  OF  NITRATE 
LEACHING  LOSSES 


M.  Mehran,  K.K.  Tanji  and  I.K.  Iskandar 


INTRODUCTION 

Land  treatment  of  municipal  wastewater  has 
gained  much  attention  in  recent  years,  primarily  be¬ 
cause  of  the  scarcity  of  water  in  arid  and  semi-arid 
regions  (Sanai  and  fchayegan  1978),  the  high  energy 
cost  of  advanced  wastewater  treatment  (Dryden  and 
Chen  1978),  the  pollution  of  surface  waters  as  a  re¬ 
sult  of  the  direct  discharge  of  sewage  effluents,  and 
the  beneficial  use  of  sewage  constituents  for  agricul¬ 
tural  production  (Bouwer  1974,  Kardos  and  Sopper 
1974).  Moreover,  the  provisions  of  sections  201  (g) 

(5)  and  313(b)(2)  of  the  Clean  Water  Act  of  1977 
require  that  land  treatment  be  fully  studied  as  a 
feasible  wastewater  treatment  alternative  before  the 
construction  of  any  facility.  Managing  a  land  treat¬ 
ment  system  requires  more  than  a  thorough  under¬ 
standing  of  the  nature  of  the  effluents,  soil  condi¬ 
tions,  and  environmental  factors.  It  also  makes  de¬ 
mands  on  the  ability  to  predict  and  control  the  be¬ 
havior  of  the  overall  system  for  optimal  long-term 
operation. 

Optimal  operation  could  be  achieved  by  resorting 
to  experimentation.  But  this  approach  would  not 
only  be  costly  and  time-consuming  but  might  not 
provide  the  generalities  we  often  need  in  extrapolating 
the  results  to  new  and  different  conditions.  Modeling, 
on  the  other  hand,  is  considered  to  be  one  of  the  most 
important  tools  for  understanding,  predicting,  and 
controlling  the  behavior  of  nitrogen  in  natural  sys¬ 
tems,  particularly  those  manipulated  by  man. 

Realization  of  the  objectives  of  the  modeling  effort 
is  perhaps  a  key  element  in  the  initial  formulation  of  a 
model.  The  existing  knowledge  of  individual  processes 
and  mechanisms,  as  well  as  of  the  interrelationships 
among  them  and  with  the  surroundings,  form  the 
foundation  for  potential  conceptualization  of  the 


problem.  Taking  the  process  of  modeling  one  step 
further,  the  degree  of  refinement  or  sophistication  of 
the  model  depends  on  the  modeler’s  ability  to  concep¬ 
tualize  and  mathematically  describe  the  events  which 
are  thought  to  be  of  prime  importance  in  attaining  the 
objectives  in  mind. 

Modeling  approaches 

The  concepts  of  systems  analysis  have  provided  the 
essential  building  blocks  of  most,  if  not  all,  modeling 
activities  in  many  diverse  disciplines  and  have  helped 
bring  an  understanding  of  the  behavior  of  systems 
under  natural  and  artificial  stresses  (Endelman  et  al. 

1972,  Van  Dyne  1970,  Schultz  1969). 

Various  modeling  approaches  have  been  taken  to 
describe  the  biophysical  and  chemical  phenomena  of 
land  application  of  wastewater.  They  can  be  broadly 
classified  as  conceptual  and  dynamic.  In  this  report  we 
refer  to  conceptual  (or  compartmental)  models  as 
those  that  attempt  to  organize  a  complex  system  into 
subsystems,  compartments,  blocks,  pools,  or  reservoirs, 
primarily  in  order  to  simplify  different  simultaneous 
processes  and  interactions  which  take  place  in  the  sys¬ 
tem.  Each  compartment  describes  part  of  the  system 
and  the  overall  system  functions  in  terms  of  fluxes  be¬ 
tween  the  subsystems  (Gordon  1969).  Although  the 
selection  of  compartments  is  to  some  extent  arbitrary, 
it  is  important  to  draw  clear  boundaries  in  accordance 
with  the  significant  physical  phenomena  governing  the 
behavior  of  the  system.  The  relative  time  scale  and 
size  of  the  compartments  and  fluxes  are  of  primary  and 
crucial  importance  in  the  formulation  of  a  compartmental 
model.  Compartmental  models  may  not  reveal  all  the 
cause-and -effect  relationships  in  the  system  but  they 
can  serve  as  tools  for  predicting  the  impact  of  long-term 
management  alternatives,  which  in  turn  can  provide 
useful  guidelines  for  decision-makers  and  managers. 


Dynamic  models  are  more  mechanism-oriented. 

They  offer  more  insight  into  the  individual  processes 
within  the  system  as  well  as  the  interrelationships 
among  different  processes.  Dynamic  models  attempt 
to  describe  the  events  as  the  best  approximation  to 
whrt  occurs  in  nature  in  comparatively  small  temporal 
and  spatial  increments.  However,  their  large  input  data 
requirements,  the  high  cost  of  obtaining  the  input  data, 
the  high  cost  of  computer  runs,  and  the  difficulties  in¬ 
volved  in  verifying  model  results  limit  their  utility. 

Existing  models 

Compartmental  models  of  nitrogen  transformation- 
have  been  applied  on  geographical  scales  ranging  from 
the  global  level  to  small  laboratory  experimental  col¬ 
umns.  On  a  global  scale,  the  emission  of  nitrous  oxide 
into  the  atmosphere  and  its  impact  on  the  ozone  con¬ 
tent  of  the  stratosphere  (Crutzen  1972)  has  raised  ques¬ 
tions  with  regard  to  the  extensive  use  of  N  fertilizers 
(Olson  1978).  Several  investigators  have  compiled 
quantitative  estimates  of  various  compartments  and 
fluxes  contributing  to  the  global  N  cycle  (Hutchinson 
1954,  Erickson  1959,  Delwiche  1970).  One  of  the  most 
recent  estimates  of  the  global  N  cycle  was  given  by 
Soderlund  and  Svensson  (1 976).  It  was  concluded 
that  the  global  N  balance  estimate  is  not  accurate  be¬ 
cause  of  uncertainties  in  estimating  over  time  the  N 
content  in  terrestrial  and  aquatic  reservoirs. 

The  intensive  use  of  N  fertilizers  in  irrigated  agricul¬ 
ture  has  caused  regional  and  local  increases  in  the  ni¬ 
trate  content  of  groundwaters  to  levels  where  they  are 
hazardous  for  human  consumption  and  other  domes¬ 
tic  purposes  (Commoner  1970).  An  increase  in  NOJ-N 
concentration  in  groundwater  in  the  Upper  Santa  Ana 
Basin,  California,  led  to  a  study  in  which  N  pools  and 
fluxes  were  estimated  in  order  to  identify  the  sources 
and  sinks  so  that  recommendations  and  guidelines 
could  be  developed  for  controlling  and  preventing 
future  problems  in  the  area  (Ayers  and  Branson  1973). 
Calculations  of  the  mass  balance  of  N  applied  to 
356,000  acres  of  the  basin  floor  showed  that  large 
amounts  of  N  per  year  are  potentially  available  for 
leaching,  accumulation  or  uptake  by  plants  in  future 
years.  This  excessive  nitrate  loading  was  attributed 
largely  to  land  disposal  of  animal  manures,  municipal 
sewage  effluents,  and  N  fertilizers. 

Aldrich  (1972),  using  another  highly  aggregated 
pool  approach  for  the  Upper  Sangamon  River  Basin, 
Illinois,  suggested  that  N  fertilizer  is  a  significant  con¬ 
tributor  to  the  nitrate  found  in  surface  waters.  Miller 
and  Smith  (1976)  followed  the  same  procedure  as  that 
used  in  the  Santa  Ana  Basin  study  and  concluded  that 
information  on  N  inputs  and  outputs  must  be  greatly 
improved  for  more  accurate  estimates  of  N  pools  and 


fluxes  within  the  rather  complex  Southern  San  Joaquin 
Valley  Basin.  Some  of  the  most  important  unknowns 
included  denitrification,  N  in  return  flows,  and  geologi¬ 
cal  formations  of  substrata  (Miller  and  Wolfe  1978). 

Tanji  (1976)  applied  a  compartmental  hydrosalinity 
model  to  the  Panoche  Drainage  District,  California,  and 
predicted  the  surface  irrigation  return  flow  to  within 
±17%  of  the  monitored  value  for  1975.  A  further  veri¬ 
fication  of  model  outputs  was  the  agreement  between 
the  calculated  total  dissolved  solids  (TDS)  of  the  soil 
solution  (7146  mg  L*1 )  and  the  measured  TDS  of  tile 
drainage  effluent  concentration  (6710-8900  mg  L'1) 
from  a  678-ha  tile-drained  farm.  Another  model  by 
Duffy  et  al.  (1975),  designed  to  predict  NO3  content 
in  the  tile  effluent,  predicted  the  tile  flow  and  the 
height  of  the  water  table  on  a  62 -ha  tile-drained  field. 

Compartmental  models  have  also  been  applied  to 
describe  the  major  changes  in  soil  N  as  a  result  of  in¬ 
coming,  outgoing,  and  internal  processes  in  greenhouse 
experiments  (Tyler  and  Broadbent  1958),  lysimeters 
(Allison  1955,  Owens  1960)  and  field  experiments 
(Boawn  et  al.  1960).  An  extensive  review  of  N  mass 
balance  studies  in  greenhouse  experiments  and  lysi¬ 
meters  was  made  by  Allison  (1965). 

Fried  et  al.  (1976)  introduced  a  long-term  concept 
for  evaluating  leaching  losses  as  a  result  of  agricultural 
practices.  It  was  based  on  the  assumption  that  the 
amounts  of  N  fertilizer  added  at  the  usual  application 
rates  in  highly  productive  soils  over  50  to  100  years  are 
so  large  relative  to  the  amounts  of  native  soil  N  that  in 
the  long  run  the  contribution  of  the  original  soil  N  to 
N  in  the  drainage  water  can  be  neglected  (Kohl  et  al. 
1978).  It  was  suggested  that  the  soil  system  will  reach 
a  steady  state  in  which  the  net  contribution  of  the  soil 
N  will  approach  zero  and  input  will  equal  output.  This 
steady-state  concept  was  then  used  to  evaluate  the  effi¬ 
ciency  of  added  N  fertilizers  and  leaching  losses. 

Based  on  the  above  concept  (Fried  et  al.  1 976)  and 
the  principle  of  mass  balance,  a  compartmental  steady 
state  model  was  formulated  to  estimate  N  emission 
from  cropped  land  (Tanji  et  al.  1977).  The  main  as¬ 
sumption  was  that  the  net  change  in  storage  of  soil  N 
is  zero  on  an  annual  basis.  The  state  variables  in  the 
model  consisted  of  nitrogen  and  water.  This  model 
was  applied  to  two  sets  of  experimental  field  data  for 
corn  monoculture.  The  model  was  applied  reasonably 
well  to  Davis,  California,  com  plots  at  the  Kearney  Field 
Station  where  the  soil  was  coarse -textured  and  low  in  re¬ 
sidual  organic  matter.  In  contrast,  the  model  grossly 
underestimated  N  leaching  losses  in  the  Davis  corn 
plots  because  both  the  contributions  from  mineraliza¬ 
tion  of  soil  organic  N  and  residual  inorganic  N  were 
neglected  in  the  model. 

This  steady-state  model  was  later  extended  so  that 
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it  would  be  applicable  to  transient  conditions,  including 
considerations  for  initial  inorganic  soil  N  as  well  as  soil 
organic  N  (Tanji  et  al.  1979).  It  was  applied  to  three 
irrigation  regimes  and  four  fertilizer  treatments.  The 
agreement  between  calculated  and  measured  values  of 
plant  uptake,  residual  inorganic  N  and  average  annual 
soil  solution  N  concentrations  was  generally  good. 

This  model  will  be  discussed  in  more  detail  later  in 
this  report. 


UNDERLYING  PRINCIPLES  FOR 
COMPARTMENTAL  MODELING 

Definition  of  the  boundaries  of  the  system  to  be 
modeled  is  perhaps  the  first  step  in  the  formulation  of 
a  compartmental  model  (National  Research  Council 
1978).  This  enables  the  modeler  to  locate  the  points 
of  inflow  to  and  outflow  from  the  system.  Identifica¬ 
tion  of  inflows  and  outflows,  as  well  as  flows  intercon¬ 
necting  the  compartments  within  the  system,  and  the 
recognition  of  major  sources  and  sinks  are  essential  com¬ 
ponents  of  model  development.  To  proceed  with  a  quan¬ 
titative  analysis,  we  have  to  define  the  size  of  the  com¬ 
partments  (subsystems),  rates  of  fluxes  across  the  com¬ 
partment  boundaries,  and  rates  of  source  and  sink  pro¬ 
cesses  contributing  to  the  state  variable  or  variables 
under  consideration.  From  the  known  values  of  fluxes 
and  compartments,  an  unknown  flux  or  compartment 
value  can  be  computed.  Depending  upon  the  nature  of 
the  input  data,  the  frequency  of  monitored  or  meas¬ 
ured  values,  and  the  rates  of  changes  in  the  sources  and 
sinks,  an  initial  time  increment  must  be  selected  which 
may  be  subject  to  change  as  the  processes  of  model  de¬ 
velopment,  calibration,  and  verification  proceed. 

Mass  balance 

One  of  the  most  fundamental  laws  of  transport 
phenomena  is  the  law  of  conservation  of  mass.  The 
concept  of  material  balance  is  not  only  a  necessary 
condition  for  formulation  of  compartmental  models 
but  for  any  model  comprising  transport  and  transforma¬ 
tion  processes.  Even  in  the  most  sophisticated  dynamic 
models,  a  test  of  the  analytical  or  numerical  schemes 
can  be  performed  by  mass  balance,  particularly  in  the 
case  of  nonlinear  problems.  Mass  balance  defines  all 
inputs,  outputs,  and  changes  in  storage  of  a  material 
within  a  bounded  region  for  a  particular  time  interval. 

In  one  space  dimension,  this  can  be  shown  mathemati¬ 
cally  in  the  following  manner: 


where  m  -  mass 

x  =  space  coordinate 
t  =  time 

s  =  sources  and/or  sinks. 

In  eq  1,  the  term  on  the  left-hand  side  represents 
changes  in  the  storage  from  time  0  to  time  t  for  the  en¬ 
tire  domain.  The  first  term  on  the  right  hand  side  de¬ 
fines  the  net  flow  of  mass  in  and  out  of  the  system  with 
boundaries  located  at  0  and  L  while  the  second  term  ac¬ 
counts  for  all  the  sources  and  sinks  contributing  to  a 
change  in  mass  m  within  the  bounded  region  during 
time  interval  t. 

Steady  state  vs  transient  state 

Compartmental  models  may  describe  systems  having 
either  steady  or  transient  behavior.  The  assumption  of 
steady  state  will  reflect  long-term  trends  and  zero  net 
change  in  storage.  In  systems  where  it  is  believed  that  a 
steady  state  prevails,  the  left-hand  side  of  eq  1  approaches 
zero  as  the  system  moves  toward  a  steady  state.  A  true 
steady  state,  which  is  rare  in  nature,  implies  that  the 
state  variables  are  not  changing  with  time.  However,  by 
selecting  a  proper  time  interval,  pseudo-steady -state  con¬ 
ditions  may  be  treated  as  steady  state  (National  Research 
Council  1978).  Kohl  et  al.  (1978)  believed  that  it  is  ad¬ 
vantageous  to  impose  the  steady-state  assumption  when 
the  compartments  are  small  compared  to  fluxes  across 
them.  One  advantage  of  imposing  a  steady-state  condi¬ 
tion  on  the  system  is  that  it  allows  us  to  compute  the 
unknown  outgoing  fluxes  from  the  incoming  fluxes  by 
equating  inputs  to  outputs.  Models  based  on  the  steady- 
state  assumption  are  generally  in  the  form  of  a  system 
of  algebraic  equations  which  could  be  solved  either  si¬ 
multaneously  or  in  sequences  (Tanji  et  al.  1977). 

Natural  systems,  particularly  for  the  time  intervals 
of  usual  interest,  are  under  transient  state  conditions. 

This  greatly  complicates  the  formulation  of  the  model 
because  of  the  nature  of  flow  rates,  space  and  time 
dimensions,  and  functional  relationships  governing  the 
transport  and  transformation  mechanisms.  For  instance, 
in  cases  where  the  rate  of  change  of  a  compartment  is 
proportional  to  the  size  of  the  compartment,  a  system 
of  simultaneous  differential  equations  can  be  derived 
which,  depending  upon  their  complexity,  can  be  solved 
analytically  or  numerically  (Mehran  and  Tanji  1974, 

Van  Dyne  1 970).  It  should  be  pointed  out  that  not  all 
transient  compartmental  models  are  formulated  in 
terms  of  differential  equations  and  it  is  possible  to  have 
a  system  of  algebraic  equations  describing  the  rate  of 
change  In  state  variables  with  time  (Tanji  et  al.  1979), 

Compartment  size  and  level  of  aggregation 

The  most  important  consideration  in  setting  up 
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compartments  is  selecting  their  size  and  the  rates  of 
fluxes  across  their  boundaries  so  that  they  can  be  meas¬ 
ured  with  a  reasonable  degree  of  accuracy  (National 
Research  Council  1978).  This  emphasizes  the  impor¬ 
tance  of  proper  selection  of  subsystem  boundaries  for 
each  compartment.  The  degree  or  level  of  aggregation 
in  a  model  depends  primarily  on  the  size  of  the  system 
and  the  objectives  of  the  modeling  effort.  Obviously, 
as  the  level  of  aggregation  increases,  the  precision  of 
the  model  outputs  decreases.  Models  should,  therefore, 
be  structured  to  have  small  compartments  relative  to 
the  fluxes  across  the  boundaries  so  that  the  changes  in 
the  size  of  the  compartments  can  be  measured.  This  is 
particularly  important  where  steady  state  conditions 
prevail,  in  which  case  the  assumption  of  constant  size 
of  the  compartment  becomes  experimentally  difficult 
to  verify.  It  is  perhaps  more  reasonable  in  these  cases 
to  resort  to  an  accurate  measurement  of  fluxes  rather 
than  compartments  (Kohl  et  al.  1 978). 

Most  conceptual  models  which  attempt  to  simulate 
the  behavior  of  N  in  soil-water-plant  systems  are  of 
the  lumped-parameter  type,  meaning  that  input-output 
relations  are  examined  only  in  terms  of  time  and  not  in 
discrete  spatial  units.  This  implies  that  the  state  vari¬ 
able  is  averaged  over  the  total  space  of  any  individual 
compartment  or  the  system. 

On  the  basis  of  available  input  data  and  existing 
data  for  validation,  a  choice  must  be  made  with  regard 
to  the  type  of  model  and  its  level  of  aggregation  to  pro¬ 
vide  the  best  methodology  for  evaluating  the  state  vari¬ 
ables  at  minimum  cost.  This  emphasizes  the  importance 
of  not  only  the  cost  of  model  development  and  valida¬ 
tion,  but  also  the  cost  of  computer  runs  for  new  condi¬ 
tions.  A  sophisticated  but  inefficient  computer  model 
has  limited  utility. 

CONCEPTUALIZATION  OF  N  TRANSPORT 
AND  TRANSFORMATION  PROCESSES 

The  main  objective  of  the  compartmental  model 
presented  here  is  to  compute  the  concentration  and 
mass  of  nitrate  leaching  losses  below  the  root  zone. 

This  general  model  is  applicable  to  chemical  fertilizer, 
municipal  wastewater,  and  animal  manure  N  inputs. 

Since  water  is  the  primary  carrier  of  mobile  N  spe¬ 
cies  and  a  medium  for  almost  all  transformation  me¬ 
chanisms,  a  knowledge  of  the  mass  balance  of  water  is 
essential  in  the  computation  for  NOJ-N  concentration 
and  mass.  The  model  is  composed  of  a  water  submodel 
and  a  N  submodel  and  is  derived  from  that  of  Tanji  et 
al.  (1979). 


Water  flow  submodel 

The  processes  contributing  to  losses  and  gains  of 
water  in  the  soil-water-plant  system  include  irrigation, 
precipitation,  evapotranspiration,  surface  return  flow, 
and  deep  percolation.  A  schematic  diagram  of  the 
fluxes  and  pools  of  the  water  flow  submodel  is  shown 
in  Figure  1 .  Knowing  the  amount  of  irrigation,  precipi 
tation,  and  evapotranspiration,  the  following  relations 
define  the  unknowns  in  terms  of  known  quantities: 


<?sw=*iae<?hv  +  0-‘)<?p  (2) 

QqX  -  C?sw  "  ^et  (3) 

<4) 

^dp^kiaeOiw+0^)Op-Oet|  (5) 

<?*lrf  =  0^iaeOiw)+flOp  +  0-0(?M  (6) 

where  Qsw  =  soil  water 


eiae  =  irrigation  application  efficiency 
<5 jw  =  irrigation  water 
a  =  runoff  coefficient 
Op  =  precipitation 
Qtt  =  leaching  loss 
0el  =  evapotranspiration 
0dp  =  deep  percolation 


Figure  I .  Schematic  diagram  of  water  flow  submodel. 


f-  deep  percolation  coefficient 

Qsirf  =  surface  irrigation  return  flow. 

All  water  fluxes  are  in  units  of  centimeters  per  year. 
Besides  the  known  quantities  mentioned  above,  three 
coefficients  must  be  determined  or  estimated  before 
the  above  equations  can  be  solved,  namely  irrigation 
application  efficiency  eiae,  precipitation  runoff  coef¬ 
ficients,  and  deep  percolation  coefficient  f. 

In  the  above  formulation,  it  is  assumed  that  the 
amount  of  evapotranspiration  is  measured  or  calculated 
so  that  leaching  losses  can  be  estimated.  This  is  the 
usual  case  for  almost  all  field  conditions  where  leaching 
losses  are  estimated  by  difference.  But  in  more  con¬ 
trolled  environments,  such  as  those  with  lysimeters 
where  drainage  is  measured,  the  above  relationships 
will  yield  a  value  for  (?et  by  a  simple  mass  balance  for 
water. 

Nitrogen  flow  submodel 

The  major  pools  and  fluxes  that  are  believed  to  con¬ 
tribute  to  N  leaching  or  accumulation  in  the  soil-water- 
plant  system  are  shown  in  Figure  2.  Applied  N  is  com¬ 
posed  of  inorganic  (/Vaif)  and  organic  (/Vaof)  fertilizers, 
N  fixation  (A/fn),  and  addition  of  N  as  a  result  of  irriga¬ 
tion  or  wastewater  (/Vjw)  and  precipitation  (/Vp).  Soil 
N  consists  of  two  measurable  pools,  namely  the  inor¬ 
ganic  N  and  organic  N  pools,  that  may  serve  as  sources, 
sinks,  or  both.  The  outputs  of  N  include  deep  percola¬ 
tion  (Wdp),  plant  uptake  (/Vhc),  denitrification  (A/d), 
and  N  in  the  surface  return  flow  (yVsjrf).  All  N  fluxes 
are  in  units  of  kg  ha'1  yr’1 . 


The  total  input  flux  of  N  to  the  soil  inorganic  pool 
can  be  described  by  the  following  equation: 

"in  =^isn+NMCFkiae-Ciw  •  <?iw  +(1  -*)Cp<?p] 

+  (1  -9)  Wajf +  (1  -d)  (1  -/"’’V, of 

+  (1-e'/Cm2,)A'ft,  +  (1-^/Cm3,)A'osn  (7) 

where  NMCF  =  N  mass  conversion  factor 

Cjw  =  N  concentration  in  irrigation  water 
Cp  *  N  concentration  in  precipitation 
g  =  gas  loss  coefficient  for  applied  inor¬ 
ganic  N 

d  =  gas  loss  coefficient  for  organic  N 
/Cml,  /Cm2,  /Cm 3=  mineralization  rate  constant 
t  -  time. 

The  inorganic  N  pool  (/Vjrz)  is  then  subject  to  denitrifi¬ 
cation  and  uptake  according  to  the  following  relations 
from  which  the  amount  of  denitrified  N  (A/d )  and  the 
amount  of  plant  uptake  of  N  (/Vhc)  can  be  obtained: 

Na=C{NJ  (8) 

and 

Nhc=h(UC)Nlrz  (9) 

where  C  is  the  denitrification  coefficient  and  h  the  crop 
uptake  coefficient.  From  the  above  relations,  the  amoun 
of  inorganic  N  subject  to  leaching  can  be  computed  as 
follows: 


Figure  2.  Schematic  diagram  of  nitrogen  flow  submodel. 
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/Vss=£(/Vir2-/Vd-/Vhc) 


(10) 


CASE  STUDIES 


where  8  is  the  leaching  fraction  defined  as 

8= - 2*L - -  (11) 


The  compartmental  model  described  above  was  ap¬ 
plied  to  two  sets  of  experimental  data  obtained  from 
1)  controlled  test  cells  located  at  CRREL,  and  2)  field 
experimental  plots  at  the  University  of  California,  Davis. 


and  thus  the  concentration  can  be  calculated  as 


N< 


fifi 


NMCF*  <?cc 


(12) 


In  cases  where  intercepted  subsurface  drainage  is  zero, 
the  mass  of  N  in  deep  percolation  /Vdp  will  be  equal  to 
as  calculated  from  eq  10,  but  otherwise  the  follow¬ 
ing  equation  should  be  used: 


*dP  =  ™M.  03) 

The  concentration  and  mass  of  N  in  surface  return  flow 
is  calculated  as 


/Vsirf  =  NMCF  [(1  -eiae)Cjw(?iw+<7Cp(?p 

+  (1-0Cse(?M]  +A/iwronp  +  /Vpronp  (14) 

where  /Vjwronp  and  yVpronp  are  the  mass  of  N  pickup  by 
irrigation  and  precipitation  runoffs,  respectively,  and 

csirf= - 0*L - .  (15) 

s,rf  NMCF  •  Qsjrf 

Residual  inorganic  /Vrjsn  and  organic  N  (/Vr0sn)  at 
the  end  of  time  t  are  calculated  as  follows: 


"risna(1-«)(/Vir,-/V/Vhc)  (16) 


+  /Vfne**m2'  (17) 

The  parameters  and  coefficients  used  in  the  model 
will  vary  according  to  the  particular  conditions  under 
study.  The  two  case  studies  given  in  the  next  section 
will  provide  more  insight  into  their  nature,  magnitude, 
and  variation.  Although  Figure  2  shows  that  there  are 
only  two  pools  of  soil  N  (Njsn  and  /Vosn),  it  should  be 
realized  that  they  are  in  turn  subdivided  into  smaller 
compartments  for  computational  purposes.  Therefore, 
the  actual  level  of  aggregation  in  the  model  is  not  as 
high  as  is  indicated  in  Figure  2. 


Case  1 -CRREL  test  cells 

In  1973  six  outdoor  cells  were  constructed  of  rein¬ 
forced  concrete  with  a  square  cross  section,  8.5  m  on 
each  side  and  1 .5  m  deep.  Cells  1 , 2  and  3  were  filled 
with  Windsor  sandy  loam  soil  and  cells  4,  5  and  6  with 
Charlton  silt  loam.  The  average  bulk  densities  of  the 
two  soils  were  1 .55  and  1 .44  g  cm-3 ,  respectively.  Sew¬ 
age  effluent  was  obtained  from  a  nearby  housing  com¬ 
munity,  given  conventional  primary  or  secondary  treat¬ 
ment,  disinfected  with  ozone,  and  applied  by  spray  irri¬ 
gation  to  the  test  cells.  Cells  1  and  6  received  5  cm  of 
secondary  effluent  per  week  during  the  growth  season 
for  five  years.  Other  cells  received  varying  amounts  of 
either  primary  or  secondary  treated  sewage  effluents 
ranging  from  5  to  15  cm  throughout  the  5-year  period 
(Iskandar  et  al.  1976,  Jenkins  et  al.  1978).  Nitrogen 
applied  to  the  test  cells  with  treated  wastewaters  was 
primarily  in  the  ammonium  form.  The  average  total  N 
concentrations  in  the  primary  and  secondary  wastewater 
used  were  26.0  and  26.9  mg  N  L'1 ,  respectively. 

Required  input  data  and  coefficients 

Total  water  inflow  to  the  cells  consisted  of  precipi¬ 
tation  and  applied  sewage  wastewater.  Water  percolat¬ 
ing  to  the  1 .5-m  depth  was  collected  and  measured, 
allowing  estimates  for  evapotranspiration  by  difference. 
Knowing  the  total  amount  of  applied  sewage  and  its  con¬ 
centration  of  N,  total  annual  N  input  was  calculated. 

Soil  organic  N  (Table  1)  was  obtained  from  soil  analysis 
data  (Iskandar  et  al.  1979).  The  amount  of  N  percolate 
for  the  period  1973-74  was  estimated  using  available 
data  on  the  amount  of  water  applied  and  the  percentage 
percolated  for  all  cells  from  1974  to  1978  (Fig.  3).  The 
time  increment  used  in  the  model  was  one  year,  begin¬ 
ning  in  June  of  each  year.  Initial  organic  N  before  the 
start  of  the  experiment  was  assumed  to  be  zero.  The 
validity  of  this  assumption  will  be  discussed  in  the  next 
section.  Since  there  was  no  surface  runoff  in  the  cells, 
the  application  efficiency  and  deep  percolation  coeffi¬ 
cients  were  considered  to  be  unity  while  the  precipita¬ 
tion  runoff  coefficient  was  taken  to  be  zero.  The  con¬ 
centration  of  N  in  precipitation  was  assumed  to  be  1 .2 
mg  L'1 .  The  amounts  of  N  fixation  and  net  mineraliza¬ 
tion-immobilization  were  also  assumed  to  be  zero. 

This  is  because  the  amount  of  applied  N  is  much 
greater  than  that  produced  by  N  fixation  and  experi¬ 
ence  has  shown  that,  in  slow  rate  systems,  immobile 
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Table  1 .  Summary  of  water  and  N  input  data  for  six  test  cells  from  1973 
to  1978  (Jenkins  etal.  1981). 


Irrigation  Applied  Organic 

water  Precipitation  Percolation  inorganic  N  soil  N 


Period 

(cm) 

(cm) 

(cm) 

(kg  ha'1 ) 

(TO3  kgha'}) 

Cell  7 

1973-74 

67 

93 

91* 

353 

44 

74-75 

149 

74 

148 

673 

62 

75-76 

92 

116 

135 

324 

59 

76-77 

93 

99 

123 

399 

71 

77-78 

67 

86 

113 

348 

71 1 

Cell  2 

1973-74 

167 

93 

210* 

928 

41 

74-75 

451 

74 

447 

2055 

58 

75-76 

216 

116 

256 

784 

70 

76-77 

87 

99 

139 

430 

72 

77-78 

41 

86 

56 

232 

72t 

Cell  3 

1973-74 

70 

93 

96* 

393 

44 

74-75 

224 

74 

260 

1072 

72 

75-76 

92 

116 

141 

361 

63 

76-77 

85 

99 

132 

429 

77 

77-78 

41 

87 

57 

233 

77f 

Cell  4 

1973-74 

71 

93 

96* 

393 

87 

74-75 

216 

74 

248 

1003 

117 

75-76 

89 

116 

165 

343 

125 

76-77 

87 

99 

126 

436 

135 

77-78 

46 

87 

63 

251 

1 35  f 

Cell  5 

1973-74 

141 

93 

183* 

602 

49 

74-75 

206 

74 

226 

956 

109 

75-76 

95 

116 

162 

358 

108 

76-77 

105 

99 

133 

533 

118 

77-78 

43 

87 

54 

231 

1 1 8 1 

Cell  6 

1973-74 

67 

93 

91* 

334 

79 

74-75 

147 

74 

177 

628 

123 

75-76 

121 

116 

179 

411 

127 

76-77 

94 

99 

122 

398 

110 

77-78 

63 

87 

97 

328 

not 

*  Estimated  from  Figure  3. 

f  Not  measured  in  1978;  value  of  previous  year  was  assumed. 


zation  equals  mineralization  over  a  short  period  of 
time  (a  few  years). 

The  overall  mass  balance  of  N  for  all  cells  after  five 
years  showed  approximately  5%  of  the  N  unaccounted 
for,  and  thus  the  denitrification  coefficient  for  the 
model  was  assumed  to  be  0.05  (the  denitrification  also 
includes  gaseous  losses  of  N).  Since  there  were  no  in¬ 
dependent  measured  values  for  the  rate  of  N  uptake  by 
grass,  the  annual  uptake  y  was  plotted  against  total 
applied  N  only  for  those  cases  which  showed  a  recovery 
of  applied  N  better  than  85%  (Fig.  4),  For  the  above  re¬ 
lationship,  the  following  expression  was  used: 


•Vmax  (*-*m)  (18) 

k  +  (x-*J 

where  ymax ,  xm,  and  k  arc  constants  with  approximate 
values  of  600,  1 50  and  400  respectively  for  the  curve 
shown  in  Figure  4. 

Using  eq  18,  the  initial  computations  of  the  model 
resulted  in  higher  annual  uptake  of  N,  particularly  for 
the  1976-76  and  1976-77  periods,  as  compared  with 
measured  values.  Because  the  experimental  data  showed 
that  potassium  (K+)  was  deficient  for  forage  growth 
during  the  above  periods  (Palazzo  and  Jenkins  1979), 

N  uptake  (Nhca)  was  adjusted  according  to  the  follow¬ 
ing  empirical  equation: 
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Totol  Applied  Water  (cm) 


Figure  3.  Percolated  water  (%)  as  a  function  of  annua!  applied  water 
for  test  cells  1  to  6  from  1974-7977. 


Figure  4.  The  relationship  between  applied  inorganic  N  and  plant  up¬ 
take  for  treatments  with  less  than  15%  unaccounted  N. 


09) 

where  0-9  is  the  optimum  value  of  the  ratio  K/N  (R ) 
in  the  plants  at  which  N  uptake  is  neither  suppressed 
nor  enhanced  by  the  amount  of  K  (National  Research 
Council  1970,  Potash  Institute  of  America  1973). 
Table  2  shows  the  ratio  K/N  in  the  plants,  the  ex¬ 
changeable  K  in  the  soil  (Iskandar  et  al.  1970,  Palazzo 


Table  2.  Measured  K/N  ratio,  soil  exchangeable  K+ 
and  R  values  for  the  model. 


Period 

K/N  in  forage 

Exch  K+  in  soil,  meq  (JOOg'1 ) 

R 

Windsor 

Chariton 

Windsor 

Chariton 

1973-74 

0.83 

0.83 

0.12 

0.09 

0.9 

74-75 

0.83 

0.73 

0.14 

0.22 

0.9 

75-76 

0.78 

0.70 

0.06 

0.07 

0.6 

76-77 

0.72 

0.80 

0.07 

— 

0.6 

77-78 

1.13 

1.58 

0.35 

0.44 

0.9 

and  Jenkins  1979),  and  assumed  values  of  R  for  the 
model.  The  increase  in  exchangeable  K  in  the  soil  and 
K/N  in  the  plants  in  1977-78  shown  in  Table  2  was 
due  to  the  application  of  300  kg  K  ha’1  as  KCI  in  May 
1977. 

A  unique  feature  of  the  model  application  in  this 
case  study  is  that  the  model  was  run  continuously  from 
1973  to  1977  for  each  cell;  i.e.  the  residual  inorganic 
soil  N  at  the  end  of  each  year  was  used  as  the  initial 
inorganic  N  for  the  subsequent  year,  except  that  the 
initial  inorganic  N  was  assumed  to  be  zero  for  the  first 
year. 

Results  and  discussion 

Figures  5-10  illustrate  the  results  of  the  model  as 
compared  with  measured  values  of  plant  uptake,  mass 
of  leaching  losses,  and  average  annual  NOJ-N  concen¬ 
tration.  Although  a  number  of  simplified  assumptions 
were  made  in  applying  this  model,  the  agreement  be¬ 
tween  calculated  and  measured  data  is  generally  good, 
particularly  with  respect  to  plant  uptake.  The  largest 
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deviations  occurred  in  concentrations  of  N  in  the 
leachate  for  1974-75  and  in  the  mass  of  N  leaching 
losses  for  1977-78.  The  best  and  the  most  consistent 
agreement  was  obtained  in  the  case  of  cells  1  and  6. 
These  cells  received  5  cm  per  week  of  secondary  sew¬ 
age  throughout  the  five-year  period  while  the  amount 
and  type  of  effluent  applied  to  the  other  cells  varied 
each  year,  which  could  be  a  cause  for  some  of  the 
deviations  between  the  measured  and  computed  re¬ 
sults.  The  good  agreement  in  most  cells  between  the 
computed  and  measured  data  for  1973-74  indicates 
that  the  assumption  of  0  kg  N  ha'1  for  initial  organic 
N  appears  to  be  valid. 

Considering  the  percentage  of  water  percolated 
(Fig.  3)  and  the  variations  in  N  uptake  with  the 
amount  of  applied  N  (Fig.  4),  it  appears  that  the  best 
use  of  N  by  grass  is  achieved  where  applied  N  and 
water  do  not  exceed  approximately  400  kg  ha'1  yr‘l 
and  170  cm  yr'1  respectively.  This  equals  approxi¬ 
mately  5  cm  per  week  of  wastewater  during  a  4-month 
growing  season  (or  application  period)  to  either  soil. 
This  resulted  in  an  annual  average  of  less  than  10  ppm 
N03-N.  The  same  conclusion  is  arrived  at  by  the 
model,  indicating  not  only  the  importance  of  applied 
water  and  N  in  the  effluent  but  also  the  amount  of 
K+  present  in  the  soil.  Reduction  of  the  K/N  ratio 
below  0.8  in  the  plants  causes  a  drastic  reduction  in 
N  uptake,  resulting  in  an  increase  of  available  N  for 
leaching.  This  is  clearly  shown  in  the  measured  leach¬ 
ing  losses  of  1975-76  and  1976-77  as  well  as  in  the 
model  results. 

Case  2-Davis  field  plots 

As  part  of  a  University  of  California  project  en¬ 
titled  “Nitrate  in  Effluents  from  Irrigated  Agriculture” 
supported  by  NSF-RANN,*  a  corn  field  trial  was  es¬ 
tablished  in  1973  at  the  Davis  campus  of  the  Univer¬ 
sity  of  California.  Three  irrigation  regimes  were  used: 
20,  60  or  100  cm  of  water  was  applied  by  sprinkler 
irrigation  during  the  growing  season  (May  through 
September),  which  corresponded  to  V3, 3/3  and  V3 
respectively  of  the  average  corn  evapotranspiration 
(ET)  for  the  1970-73  period  (Biggar  et  al.  1974).  All 
plots  were  pre-irrigated  prior  to  planting  so  that  the 
surface  3  m  of  the  soil  contained  100  cm  of  water. 
Fertilizer  was  applied  at  four  rates:  0, 90, 1 80  and 
360  kg  N  ha'1  of  (NH4)2S04  at  planting  time  at  a 
5-cm  depth.  The  plots  were  instrumented  with  por¬ 
ous  ceramic  cups  to  extract  soil  solutions  at  depths 
of  0.3, 0.6, 1 .2, 1 .8,  2.4  and  3.0  m.  After  harvest, 
soil  cores  were  taken  at  a  3-m  depth  to  determine 
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total  inorganic  and  organic  N  present.  The  details  of 
the  experimental  layout  and  procedures  were  reported 
by  Biggar  et  al.  (1974).  This  model  was  applied  to  two 
years  of  data,  1974-75  and  1975-76  (Tanji  et  al.  1979). 

Required  model  input  data  and  coefficients 

Annual  fluxes  of  water  used  as  input  to  the  model 
are  summarized  in  Table  3. 

The  amounts  of  inorganic  and  organic  N  obtained  in 
October  of  1974  and  1975  after  crop  harvest  are  taken 
as  initial  values  for  the  1974-75  and  1975-76  periods. 
The  data  for  all  treatments  are  given  in  Table  4. 

The  coefficients  required  by  the  water  submodel 
and  the  N  submodel  are  given  in  Table  5.  Since  the 
water  percolating  past  the  3-m  depth  was  not  directly 
measured  as  it  was  in  the  CRREL  test  cells,  evapotrans¬ 
piration  loss  C?el  was  computed  by 

Q,i^(DpJ^b(Dptl)  (20) 

where  e  and  b  are  the  corn  and  fallow  ET  coefficients, 
respectively.  Dpet  and  0pel  are  potential  ET  for  grass 
for  the  growing  and  nongrowing  seasons,  respectively, 
and  are  obtained  from  the  following  relations: 


Opet  =  Kp^pan) 

^pet  “  ^p(^pan) 


(21) 


(22) 


where  Kp  is  potential  ET  coefficient  for  grass  and  Epan 
and  £'pan  are  pan  evaporations  for  the  growing  season 
and  the  fallow  period,  respectively.  Rainfall  data,  pan 
evaporation  data  and  ET  coefficients  were  obtained 
from  a  nearby  agroclimatic  station  and  weighing  ly- 
simeter.  Using  eq  20,  the  annual  amount  of  leaching 
past  the  3-m  soil  depth  was  computed  from  eq  3. 

The  denitrification  coefficient  C  was  obtained  by 
sensitivity  analysis  of  the  %  ET  irrigation  regime  re¬ 
ceiving  360  kg  ha'1  of  inorganic  N  (Tanji  et  al.  1979). 
The  uptake  coefficients,  obtained  by  Broadbent  (1976) 
from  isotopically  labeled  N  data,  were  0.58,  0.56,  and 
0.39,  respectively,  for  90, 180  and  360  kg  ha'1  of  ap¬ 
plied  N.  An  extrapolated  value  of  0.6  for  the  uptake 
coefficient  for  the  0-kg  N  ha’1  treatment  overestimated 
the  amount  of  uptake  (Tanji  et  at.  1979),  By  a  trial 
and  error  method,  it  was  found  that  a  value  of  0.35 
for  the  uptake  coefficient  produced  the  best  results. 
From  the  observed  N  deficiency  of  the  corn  in  the 
0-kg  N  ha'1  treatments,  it  appears  that  there  is  a 
threshold  value  of  10  mg  L'1  NOJ-N  concentration  in 
the  soil  solution  below  which  the  uptake  efficiency 
decreases  sharply. 
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Table  3.  Summary  of  water  input  data  (cm  of  H*0)  for  the 
three  irrigation  treatments  in  the  1975  (Oct  1974-Sept  1975) 
and  1976  (Oct  1975-Sept  1976)  Davis  com  plots. 


1974-1975 

1975-1976 

Item 

7,  er 

7,  st 

7,  ST 

7,  ST 

*/,  ET 

Irrigation  water 

18.9 

59.3 

100.7 

20.2 

60.5 

102.7 

Preirrigation* 

17.2 

5.4 

6.2 

40.5 

28.3 

5.9 

Precipitation 

46.5 

46.5 

46.5 

17.5 

17.5 

17.5 

Pan  evaporation 
(growing  season) 

110,4 

110.4 

110.4 

100.8 

100.8 

100.8 

Pan  evaporation 
(fallow  period) 

76.4 

76.4 

76.4 

99.5 

99.5 

99,5 

♦Preirrigation  prior  to  planting  the  1976  corn  crop  was  comparatively 
high  due  to  a  dry  winter. 


Table  4.  Initial  inorganic  and  organic  soil  N  (kg  ha-1  in  3-m  profile)  for  four 
fertilizer  and  three  irrigation  treatments  in  the  Davis  com  field  plots. 

From  analyses  of  soil  N  after  harvest  in  October  1974  and  October  1975. 


Period 

Nitrogen 
applic.  rate 

7» 

ET 

V. 

ET 

V, 

ET 

(kg  ha'1 ) 

fnorg.  N 

Org.  N* 

fnorg.  N 

Org.  N* 

inorg.  N 

Org.  N’ 

1974-1975 

0 

122 

22,200 

123 

22,200 

121 

21,811 

89.6 

132 

22,200 

104 

22,200 

119 

22,200 

179.2 

198 

22,200 

153 

22,200 

153 

22,200 

358.4 

364 

22,200 

331 

22,182 

342 

21,711 

1975-1976 

0 

110 

21,947 

120 

21,936 

117 

22,454 

89.6 

133 

22,438 

138 

21,204 

137 

20,991 

179.2 

184 

22,251 

139 

21,052 

134 

21,824 

358.4 

423 

23,081 

409 

295 

22,182 

♦Where  organic  soil  N  was  not  determined,  an  average  value  of  22,200  kg  N  ha'1  was 
assumed. 


Table  5.  Required  input  parameters  and  coefficients  for  the 
water  and  the  N  submodels. 


Item 

Symbol 

Value 

Irrigation  application  efficiency 

e\ae 

1.0 

Precipitation  runoff  coefficient 

a 

0.2 

Potential  ET  coefficient  for  grass 

kp 

0.76 

Corn  ET  coefficient 

e 

0.82 

Fallow  soil  ET  coefficient 

b 

052 

Deep  percolation  coefficient 

f 

1.0 

N  concentration  in  irrigation  water  (mg  L‘‘ ) 

ciw 

0.2 

N  concentration  in  rainwater  (mg  L'1 ) 

CP 

1.2 

Gas  loss  coefficient  for  inorganic  N 

9 

0.0 

Gas  loss  coefficient  for  organic  N 

d 

0.0 

Mineralization  rate  for  organic  N  (week’1 ) 

*m1>  *m2-  *m3  »«* 

Denitrification  coefficient 

c 

0.15 

12 
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Figure  11.  A  comparison  between  meas¬ 
ured  and  computed  values  of  N  uptake  as 
a  function  of  applied  N  for  the  5/3  ET 
treatment. 


Results  and  discussion 

Figures  11-19  show  the  comparison  between  the 
measured  and  calculated  data  for  N  uptake,  residual 
inorganic  soil  N  and  average  annual  N  concentration 
in  the  soil  solution  at  the  3-m  depth  for  the  s/3,  3/3 
and  ‘/a  ET  treatments  of  1974-75  and  1975-76.  The 
generally  good  agreement  between  the  measured  and 
calculated  data  indicates  that  the  assumed  values  of 
denitrification  and  mineralization  are  reasonable  for 


5/3  ET  Regime 


Figure  12.  A  comparison  between  meas¬ 
ured  and  computed  values  of  residua!  inor¬ 
ganic  N  as  a  function  of  applied  N  for  the 
s/3  treatment. 


the  soil  and  environmental  conditions  in  the  Davis 
plots. 

Table  6  gives  mass  emission  of  N  leaching  losses 
past  the  3-m  depth  for  the  3/3  and  sf3  treatments. 

Zero  N  emission  is  calculated  for  the  */3  ET  treatments 
because  there  was  no  downward  flow  of  water  at  the 
3-m  depth  due  to  the  evapotranspiration  losses  exceed¬ 
ing  the  effective  water  input  at  the  soil  surface. 


Table  6.  Computed  mass  emissions  of  N  (kg  N  ha’1  yr*1 )  below 
the  root  zone  (3  m)  for  the  1974-1975  and  1975-1976  corn 
field  trials  at  the  Davis  site. 


Ammonium  sulfatt 
fertilizer  rate 
(kg  ho") 

1974-1975* 

1975-1976* 

%  ET 

’/,  ET 

V,  ETt 

%ET 

%  et 

7,  ETf 

0 

28.0 

2,4 

0 

13.2 

0.52 

0 

89.6 

37.0 

3.2 

0 

20.2 

0.78 

0 

179.2 

53.1 

4.8 

0 

26.8 

t  .03 

0 

358.4 

133.4 

12.0 

0 

66.7 

2.95 

0 

*5/i  ET,  3/j  ET,  and  '/3  ET  are,  respectively,  100-,  60*,  and  20-cm  irriga¬ 
tions  during  the  growing  season. 

fN  was  not  leached  past  the  3*m  depth  because  of  upward  flow  of  water. 
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5/3  ET  Regime 


Applied  N  { kg  N  ho  *) 


Figure  13.  A  comparison  between  meas¬ 
ured  and  computed  values  of  NO 3  concen¬ 
tration  as  a  function  of  applied  N  for  the 
%  ET  treatment. 


3/3  ET  Regime 


Applied  N  (kg  N  ho*1) 

Figure  14.  A  comparison  between  meas¬ 
ured  and  computed  values  of  N  uptake  as 
a  function  of  applied  N  for  the  %  ET 
treatment. 


3/3  ET  Regime 


•l 

Applied  N  (kg  N  he  1 


Figure  15.  A  comparison  between  meas¬ 
ured  and  computed  values  of  residua!  inor¬ 
ganic  N  as  a  function  of  applied  N  for  the 
%  ET  treatment. 


3/3  ET  Regime 


Applied  N  (kg  N  ho'1) 


Figure  16.  A  comparison  between  meas¬ 
ured  and  computed  values  of  NO\  con¬ 
centration  as  a  function  of  applied  N  for 
for  the  %  ET  treatment. 
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1/3  ET  R«gim« 


o  Colcd. 


Figure  17,  A  comparison  between  meas¬ 
ured  and  computed  values  of  N  uptake  as 
a  function  of  applied  N  for  the  V3  ET 
treatment. 


1/3  ET  &•*■**• 


Figure  18.  A  comparison  between  meas¬ 
ured  and  computed  values  of  residua!  in¬ 
organic  N  as  a  function  of  applied  N  for 
the  V3  ET  treatment. 


1/3  ET  Rffltm# 


Figure  19.  A  comparison  between  meas¬ 
ured  and  computed  values  of  NO\  concen¬ 
tration  as  a  function  of  applied  N  for  the 
V3  ET  treatment. 


ROLE  OF  COMPARTMENTAL  MODELS 
FOR  PLANNERS  AND  MANAGERS 

The  precision  of  the  results  obtained  from  a  compart- 
mental  model  is  a  function  of  the  quality  of  available 
data  and  the  level  of  aggregation  of  compartments  and 
fluxes.  As  the  level  of  aggregation  increases,  the  preci¬ 
sion  decreases  (National  Research  Council  1978),  and 
of  course  the  choice  of  that  level  depends  on  the  kinds 
of  questions  that  are  being  asked  by  the  modeler.  With 
proper  selection  of  time  increment  and  size  of  pools 
and  fluxes  giving  rise  to  a  meaningful  averaging  of  events 
over  time  and  space,  compartmental  models  can  be  as 
precise  as  dynamic  simulation  models.  The  use  of  finely 
tuned  dynamic  models  for  management  purposes  may 
introduce  difficulties  as  a  result  of  1)  insufficient  infor¬ 
mation  regarding  the  processes  and  their  interrelation¬ 
ships,  2)  problems  in  obtaining  required  input  data, 

3)  high  cost  of  model  application,  and  4)  knowledge  of 
numerical  analysis  and  sophisticated  solution  techniques. 
On  the  other  hand,  compartmental  models,  according 
to  Levins  (1966),  sacrifice  realism  for  the  sake  of  gen¬ 
erality. 

Verification  of  the  compartmental  model  by  the  two 
case  studies  presented  in  this  chapter  as  well  as  those 
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models  cited  in  the  literature  (see  Existing  Models)  il¬ 
lustrate  the  potential  capabilities  of  compartmental 
models  in  predicting  the  behavior  of  N  in  soil-water- 
plant  systems  under  artificial  as  well  as  natural  stresses. 
The  uncertainties  that  exist  in  measuring  some  of  the 
pools  and  fluxes  of  N  as  well  as  the  rates  of  transforma* 
tion  impose  limitations  on  the  use  of  the  models  (Kohl 
etal.  1978,  National  Research  Council  1978).  The  fact 
that  these  limitations  should  be  evaluated  for  the  par¬ 
ticular  conditions  for  which  the  models  are  used  em¬ 
phasizes  the  importance  of  team  efforts  in  applying 
the  models.  It  is  the  cooperation  of  modelers,  re¬ 
searchers  and  decision-makers  that  sets  the  stage  for 
a  continuous  feedback  from  and  input  to  the  model 
so  that  the  model  can  be  improved.  In  other  words, 
the  model  must  be  considered  dynamic  in  nature,  as 
in  a  real  system,  to  provide  a  forum  for  future  revi¬ 
sions  and  modifications,  and  hopefully  increased 
utility. 
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APPENDIX  A:  DESCRIPTION  AND  PROGRAM  LISTING  OF 
THE  COMPARTMENTAL  NITROGEN  MODEL 


The  program  for  the  compartmental  model  consists  of  only  one  main  program.  All  input  data  are  given 
in  free  format.  The  units  for  all  N  and  water  values  are  in  kg  ha"1  yf 1  and  cm  yr_1 ,  respectively.  Input 
data  are  as  follows: 


Record  (1): 

EPAN  =  Pan  evaporation  (growing  season),  cm 
EPANP  =  Pan  evaporation  (nongrowing  season),  cm 
KP  =  potential  ET  coefficient  for  grass 
B  =  fallow  soil  ET  coefficient 

Record  (2): 

YMAX  =  maximum  uptake  of  N 
XMIN  =  minimum  applied  N  for  which  uptake  is  zero 
KM  =  constant  for  curvature  of  the  uptake  vs  applied  N 


Record  (6): 

G  =  gas  loss  coefficient  for  inorganic  N 
NAOF  =  applied  organic  N 

D  =  gas  loss  coefficient  for  organic  N 
KM  I  =  mineralization  rate  for  applied  organic  N 
NMCF  =  N  conversion  factor 

Record  (7): 

NFN  =  amount  of  N  fixed 

KM2  =  mineralization  rate  for  fixed  N 


Record  (3): 

EIAE  =  irrigation  application  efficiency 
A  =  precipitation  runoff  coefficient 
F  =  deep  percolation  coefficient 

Record  (4): 

NISN  =  resident  inorganic  soil  N 
Record  (5): 

CIW  =  N  concentration  in  irrigation  water 
CP  =  N  concentration  in  precipitation 


Record  (8): 

C  =  denitrification  coefficient 
KM3  =  mineralization  rate  for  resident  soil  organic  N 

Record  (9) : 

NCELL  =  cell  number 
NY  =  period  of  growth 
QIW  =  amount  of  applied  wastewater 
QP  =  amount  of  precipitation 
QPER  =  amount  of  water  percolated 
NAIF  =  applied  inorganic  N  at  the  soil  surface 
NOSN  =  resident  organic  soil  N 
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APPENDIX  B:  A  TRANSIENT  COMPARTMENTAL  MODEL  FOR  PREDICTION  OF 
NITRATE  LEACHING  LOSSES  AND  NITROGEN  UPTAKE  BY  GRASS  IN  SLOW 
INFILTRATION  SYSTEMS  OF  LAND  APPLICATION  OF  WASTEWATER 
C** 

C**  DEFINITION  OF  INPUT  VARIABLES  : 

C  *  *  Q  I W  =  IRRIGATION  WATER 

C**  EIAE  =  IRRIGATION  APPLICATION  EFFICIENCY 

C**  QP  r  PRECIPITATION 

C**  A  =  PRECIPITATION  RUNOFF  COEFFICIENT 
C**  EPAN  =  PAN  EVAPORATION  FOR  GROWING  SEASON 
C  *  *  EP  ASP  =  PAN  EVAPORATION  FOR  NONGROWING  SEASON 
C**  KP  =  POTENTIAL  FT  COEFFICIENT  FOR  GRASS  DURING 
C  *  *  GROWING  SEASON 

C**  KPP  =  POTENTIAL  ET  COEFFICIENT  FOR  GRASS  DURING 

C**  NONGROWING  SEASON 

C**  F  -  DEEP  PERCOLATION  COEFFICIENT 

C**  E  =  CROP  ET  COEFFICIENT 

C**  B  =■  FALLOW  ET  COEFFICIENT 

C**  N  =  NITROGEN 

C  *  *  CIW  =  N  CONCENTRATION  IN  IRRIGATION  WATER 

C**  ClWRNP  =  N  PICKUP  CONC.  BY  IRRIGATION  RUNOFF  WATER 

C**  CP  =  N  CONC.  IN  PRECIPITATION 

C**  CPRONP  =  N  PICKUP  CONC.  BY  PRECIPITATION  RUNOFF  WATER 

C**  G  =  GAS  LOSS  COEFFICIENT  FOR  INORGANIC  N  FERTILIZER 

C**  NAIF  =  APPLIED  INORGANIC  N  FERTILIZER 

C**  NAOH  =  APPLIED  ORGANIC  NITROGEN  FERTILIZER 

C**  D  =  GAS  LOSS  COEFFICIENT  FOR  ORGANIC  N  FERTILIZER 

C**  KM1  *,  KM?  H  K  M3  =  MINERALIZATION  RATE  CONSTANTS 

C**  T  =  TIME 

C**  NFN  =  FIXED  N 

C**  NOSN  =  RESIDENT  ORGANIC  SOIL  N 

C**  NISN  =  RESIDENT  INORGANIC  SOIL  N 

C**  C  =  LOSS  COEFFICIENT  FOR  DENITRIFICATION 
C**  H  =  HARVESTED  CROP  COEFFICIENT 

C** 

C** 

C  *  * 

REAL  KP»KPP*NAIF*NAOF«KMl*NFN *KM2* NOSN *KK3 *N I SN ♦ 

1NMCF  tNElWtNl W RO « NP » N EP . NPR 0 , N IFGL.NE AIF  fNOFGL« 
2NEA0F»NKA0F*NA0FM»NKFN<NFNM»NRCSN»NGSNM*NKSN« 
3NIR2»ND*NEIRZ»NHC»Nf  OR  Z«L  t NR  I SN tNL  L  »NDP  *NSDw  *NS ! R F 
CALL  ATTDEV<6*8*2*f,6) 

NR  =  5 
NW  =  6 

CALL  CONTRL<l»*INPUT*»NR> 

CALL  CONTR.L  (2»,0UTPUT*tNW) 

C**  HYDROLOGY  SUBMODEL 

C  *  * 

C**  INPUT  DATA  AND  COEFFICIENTS  FOR  HYOROLOGY  SUBMODEL 
RtAD(NR**)£IAE.A»KP»KPP*F»E*f< 

RF  AD(NR,*)EPANfLPANP 
REAQ(NR»*)GIW»QP 

C**  COMPUTATION  OF  DEEP  PERCOLATION  AND  SURFACE  RUNOFF 
QEIW  =  LIAE  *  U 1 W 
QlWRQ  =  (1.0-EIAE)*  QIW 
Q PR  0  =  A  *  OP 
QEP  =  (1.0-A)*  QP 
QSw  =  QEIW  ♦  QEP 
DPET  =  KP  *  EFAN 
DPETP  =  KPP  *  EPANP 
QET  =  £  *  DPET  ♦  3  «  DPETP 
QLL  =  QSW  -  QET 
OOP  =  F  *  QLL 
Q  SO  W  =  (1.0-F)  «  QLL 
QSIRF  =  QIwRO  ♦  Q PRO  *  OSDw 
C**  NITROGEN  SUBMODEL: 

C  *  * 
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INPUT  DATA  AND  COEFFICIENTS  FOR  NITROGEN  SUBMODEL 
R  L  A  D  !  N  P  *  *  )CI  U  ♦  C  I  W  R  N  P  tCPtCPRONP 
RE  A  D  !  N  R  »  *  )  G  *  N  A  0  F  »  D  ♦  K  M  1  »  T  «  N  M  C  F 
READ!NR,*>NFN,KM2 
RLAD!NR.*)NAIF 
REAO!NR**)NOSN 
RE AD1NR  ♦  *  >KM3 
R  L  AO  <  NR  *  *  )  N I  SN 
READC\Rt*)C,H 

COMPUTATION  OF  N  SPECIES  IN  SOIL  8,  SURFACE  W  AT  FRS 
IRRIGATION  WATER  : 

NE  I  *  =  CIW  *  QE I W  *  NMCF 
C  I « R  0  =  CIW  ♦  CIWRNP 
NIWRQ  =  CIWRO  *  3IWR0  *  NMCF 
PRECIPITATION  : 

NP  -  CP  *  QP  *  NMCF 

NEP  =  CP  *  QEP  *  NMCF 

C PR  0  =  CP  ♦  CPRONP 

NPRO  =  CPRO  *  OPRO  *  NMCF 

APPLIED  INORGANIC  N  FERTILIZER  : 

NIFGL  =  G  *  NAIF 
NEAIF  =  I1.0-G)  *  NAIF 
APPLIED  ORGANIC  N  FERTILIZER  : 

NOFGL  =  D  *  NAOF 

NEAOF  =  C1.0-D)  *  NAOF 

NRAOF  =  NEAOF  *  EXP<-KM1  *  T) 

NAOrM  =  NEAOF  -  NRAOF 
FIXED  NITROGEN  : 

NRFN  =  NFN  *  EXP!-KM2  *  T) 

NFNM  =  NFN  -  NRFN 

RESIDENT  ORGANIC  SOIL  NITROGEN  : 

NROSN  =  NOSN  *  EXP ( -KM  3  *  T) 

NOSNM  =  NOSN  -  NROSN 

RESIDUAL  ORGANIC  SOIL  N  AFTER  HARVEST  : 

NRSN  =  NROSN  ♦  NRAOF  ♦  NRFN 
NITROGEN  INFLOW  TO  ROOT  ZONE  : 

NIRZ  =  NEIW+NEP*NEAIF+NAOFM+NISN*NFNM*NOSNM 
DENITRIFICATION  AND  UNACCOUNTED  LOSSES  I 
ND  =  C  *  NIRZ 

EFFECTIVE  N  INFLOW  TO  ROOT  ZONE  : 

NEIRZ  =  NIRZ  -  ND 
HARVESTED  CROP  UPTAKE  : 

NHC  =  H  *  NEIRZ 

RESIDUAL  INORGANIC  SOIL  N  AFTER  HARVEST  Z 
NEORZ  =  NEIRZ  -  NHC 
L  =  GLL/ ( QE  IU  ♦  QEP) 

NRISN  =  tl.O-L)  *  NEORZ 
LEACHING  LOSSES  FROM  ROOT  ZONE  Z 
NLL  =  L  *  NEORZ 
CLL  =  NLL/ ( QLL  *  NMCF) 

NITROGEN  IN  RETURN  FLOWS  Z 
C  DP  =  CLL 
CSDW  =  CLL 

NDP  =  COP  *  Q DP  *  NMCF 
NSC W  =  CSDW  *  QSDW  *  NMCF 

CSIRF  =  ! C I WR  0*QI WRO+CPRO  *QPR  O+CSDW  *QSDV ) /QS 1R  F 
NSIRF  =  CSIRF  *  QSIRF  *  NMCF 


WRITE! NU *11) 

FORMAT < 1H1 « *INPUT  DATA  AND  COEFFICIENTS  FOR  HYDROLOGY  SUBMODEL*//) 
WRI TK(NWtl?> 

FORMAT (  *  QIW  LI AE  QP  A  EPAN  KP 

1  EPANP  KPP  F  E  B*//> 

WRITE!  NWt 13 )Q IW»EI AE*QP» AtEPANtKP«EPANP*KPP,F,E»B 
F0RMAT(?X*F7.2«3X.F4.2*4XtF5.1«*>X*F4.2*3X*F7.2t4X* 
lF4.2*2X*F7.2f3X«F4,2»6X*F3.1 t5XtF4.2*4X*F4.2//> 

WRITE! NW»14) 

FORMAT  C1H0** INPUT  DATA  AND  COEFFICIENTS  FOR  N  SU6M00EL*//) 

WRITE! NW«15) 


I 


r 


15  FORMAT ( 1  H  0  *  G  WAIF  NAOF  0  KM1  T  NF 

IN  KM2  NQSN  KM3  NISN  C  H  NMCF’ 

2/) 

WRITE! Ny»16)GfNAlF»NA0F*D»KMl»T»NFNfKM2tN0SN»KM3*NISN»C*H»NMCF 
16  F0RMAT!2X*F3.1»2XtFfe.2*3X«F6.2*3X*F6.2*3X»Fif.2*3X*FE.2*3X*F6.2*3X* 

lF6.4*3X*F7.1*3X*F6.5»3X*F6.2«3XtF6.2«3X<F6.2*3XtF6.2//) 
WRITE!NW*1C)CIW 

10  FORMAT ! 1H0,23X* ’N  CONC.  IN  IRRIGATION  W A TER = * « F 7 . 2 ) 

WRITE! NW»9)CP 

9  FORMAT !1H0?20X» »N  CONC.  IN  PR E C I PI T A T I ON= • ♦ F 7. 2 // > 

C” 

wR I TE( NW  » 1 7 ) 

17  FORMAT ( 1H1 *  'OUTPUT  VARIABLES  OF  HYDROLOGY  SUBMODEL*/) 

WRITE! NW*13)bEIW 

13  FORMAT 1 1h0»20X»*EFFLCTIVE  IRRIGATION  WATER  IN  CM=’fF7. 

12  ) 

WRITE!NV«19)QIwRQ 

19  FORMAT (1HO»23X»*IRRIGATION  WATER  RUNOFF  IN  CM=*,F7.2) 

WRITE! NW»20)QEP 

20  FORMAT ( 1H0  ♦  20  X  » ’EFFECT IVE  PRECIPITATION  IN  CM=*,F7.2) 

WRITE! NW»21  )QPRO 

21  FORMAT (1HC*20X**PRECIPITATI0N  RUNOFF  IN  CM=’«F7.2) 

WRITE !NW»22)GSW 

22  FORMAT (1HO»20X»’SOIL  WATFR  IN  CM=»,F7.2) 

WRITE! NW»?3)&ET 

23  FORMAT !1H0*20X»  ’EVAPOTRANSPIR  ATION  IN  CM  =  ’»F7.2) 

WRITE!NW»24  )QLL 

24  FORMAT  1 1H0«20X« ’LEACHING  LOSSES  IN  CM=’,F7.2> 

WRITE!NW*25)UDP 

25  FORMAT !1H0*20X«  ’DEEP  PERCOLATION  IN  CM=’,F7.2) 

URITE!NW»26)OSDW 

26  FORMAT!1HO,20X, ’SUBSURFACE  DRAINAGE  IN  CM=»fF7.2) 

WR  ITE <NW»27)GSIRF 

27  FORMAT  UH0»20X.  ’SURFACE  RETURN  FLOW  IN  CM=’«F7.2//) 

WRITE <NWt 30) 

30  FORMAT ! 1H0* ’OUTPUT  VARIABLES  OF  N  SUBMODEL’//) 

WRITE(NW»31) 

31  FORMAT ! 1 H  0  »  5  X  » * IRRIGATION  WATER  :*/> 

URITE(NW»32)NEIW 

32  FORMAT! 1H0*20X,’N  MASS  FROM  EFFFCTIVE  IRRIGATION  W AT ER  =  •  *F 7 . 2  ) 
URITF!Nw»33)NWR0 

33  FORMAT  I  1H0  »20X*  *N  MASS  FROM  RUNOFF  W  ATF.R  =  •  «  F  7. 2  ) 

C** 

WRITE! NW»36) 

36  FORMAT !1HO»5X**PRECIPITATION  :»/) 

WRI TE!NW»37)NEP 

37  FORMAT  <1H0*20X*’N  MASS  FROM  EFFECTIVE  PR  EC  I P I T A T I  ON- ’  »F 7 .2  ) 
WRITE(NW»38)NPR0 

38  FORMAT  1 1H0«20X»  »N  MASS  IN  PRECIPITATION  RUNOFF - •  « F 7 . 2 ) 

WRI TE ( NW  «41  ) 

41  FORMAT !  1H0 »5X t’APPLIED  INORGANIC  &  ORGANIC  FERTILIZER  I’/) 
WRITC<NW,42)NCAIF 

42  FORMAT  1 1H0»20X» ’EFFECT IVE  APPLIED  INORGANIC  NITROGEN  FERTILIZER 
1=’F7.2> 

WRI rt!NW*43)NEA0F 

43  FORMAT  1 1H0*20X»  ’EFFECT IVE  APPLIED  ORGANIC  N  FERT IL I ZER= * *F7.2 ) 
WRITE! NW*44)NA0FM 

44  FORMAT ( 1H0«20X« ’APPLIED  ORGANIC  N  FERTILIZER  M I NER AL I ZED= * *F 7 .2 > 
WRITE!NW»45)NRISN 

45  FORMAT! 1  HO » 20 X ♦ »R ES  I  DU AL  INORGANIC  N  AFTER  H AR VEST  =  • *F 7 . 2 ) 

WRITE ! NW  *46 ) 

46  FORMAT ! 1 HO ♦ 5X  •  • F I  XL D  N  J  RLSI DENT  ORGANIC  SOIL  N  JRESIOUAL  ORGANIC 
1  SOIL  H  AFTER  HARVEST’/) 

WRITE(NV»47)NFNM 

47  FORMAT ! 1H0»20X*  ’FIXED  N  M I NE R AL 1 ZE D=  •  «F 7 .2 ) 

WRITE(NW»48)NGSNM 

48  FOR  MAT !1H0*20X*’RESIDENT  ORGANIC  SOIL  N  M I NER A L  IZED=  •  « 

1F7.2) 

WRITE!NW»49)NRSN 

49  FORMAT !1HO*20X*»RESIDUAL  ORGANIC  SOIL  N  AFTER  H ARVEST= * «F8 •  1  > 
WR1TE(NW*50  ) 


r 


50  FORMAT !1H0»5X*‘DENITKIF ICATION  &  UNACCOUNTED  LOSSFS:*/) 
WRITE<NW»51)ND 

51  FORMAT !1HO*20X**DENITRIFICATION  &  OTHER  LGSS ES = • « F7. 2 ) 

WR I TE ! N W  *  53 ) 

53  FORMAT (1H0»5X»*EFFELT1VL  INFLOW  TO  ROOT  ZONE  5  HARVESTEn  CROP  UPTA 

ike:*/) 

WRITE! N  W»52)NIRZ 

52  FORMAT ! 1H0 . 20X» ‘MINER Al IZED  N  INPUT  TO  ROOT  Z0NE=»fF7.2) 

WRITEINW. 54)NEIRZ 

54  FORMAT  C 1H0*2DX*  ‘EFFECT 1VE  N  INFLOW  TO  ROOT  ZCNE=*»F7.2 
1) 

WRITE< NW«55)NHC 

55  FORMAT ( lH0t20X, ‘HARVESTED  CROP  UPT A K E  =  *  * F 7 . 2 ) 

WRITE!NW» 56) 

56  FORMAT ! 1H0 *5X  ♦ ‘LEACHI NG  LOSSES  FROM  POOT  ZONE:*/) 

WRITE!NW»57)NLL 

57  FORMAT ! 1H0 , 20X, ‘N  LEACHING  LO SSES= • « F 7  .  ? ) 

WRITE(NW»58)CLL 

58  FORMAT !1H0*20X. *N  CONC.  IN  LE  ACHATE »  MG  N/L=*»F7.2) 

WRITE«NW*6C) 

60  FOR  MAT (1H0»5X**N  IN  RETURN  FLOWS:*/) 

WRITEINW *61)NDP 

61  FORMAT ( 1  HO 1 20  X i  ‘MASS  EMISSION  THRU.  DEEP  PE R COL  AT  I ON  =  • 

1  .F7.2) 

WRITE (NW»fc2)NSDW 

62  FORMAT ( 1H0»20X»‘MASS  EMISSION  THRU.  COLLECTED  SUBSURFACE  DR  A I NAGE  = 
1  •  *  F  7 .2) 

WRITEINW *63)NSIRF 

63  FORMAT! 1H0.20X, ‘MASS  EMISSION  THRU.  SURFACE  IRRIGATION  RETURN  FLOW 
1-»*F7.P) 

WRITE! NW«64)CSIkF 

64  FORMAT ( 1H0*20X»*N  CONC.  IN  SURFACE  IRRIGATION  RETURN  FLOW=‘*F7.2) 
CALL  C0NTRL!4**INPUT‘»NR) 

CALL  C0‘JTRL!4,‘0UTPUT*.NW) 

CALL  EXIT 
END 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 


Mehran,  M. 

Evaluation  of  a  compartmental  nitrogen  model  for  pre¬ 
diction  of  nitrate  leaching  losses  /  by  M.  Mehran,  K.K. 
Tanji  and  I.K.  Iskandar.  Hanover,  N.H.:  U.S.  Cold  Regions 
Research  and  Engineering  Laboratory;  Springfield,  Va. : 
available  from  National  Technical  Information  Service, 

1981. 

v,  28  p.,  illus.;  28  cm.  (  CRREL  Report  81-23.  ) 

Prepared  for  Office  of  the  Chief  of  Engineers,  by  Corps 
of  Engineers,  U.S.  Army  Cold  Regions  Research  and  Engineer¬ 
ing  Laboratory  under  CWIS  31633. 

Bibliography:  p.  16. 
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